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ABSTRACT 

We study one (or two) matrix models modified by terms of tfie form 
g{p{P)Y + g'^p'iO))"^, wliere tlie matrix representation of the puncture 
operator P and the one of a scahng operator O are denoted by p{P) and 
p'{0) respectively. We rewrite the modified models as effective theories of 
baby universes. We find an upper bound for the gravitational dimension 
of O under which we can fine tune the coupling constants to obtain new 
critical behaviors in the continuum limit. The simultaneous tuning of 
g and g' is possible if the representations p{P) and p'{0) are chosen so 
that the non-diagonal elements of the mass matrix of the effective theory 
vanish. 
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1 Introduction 



Matrix models have played an important role in the theory of quantum gravity 
in two dimensions. The existence of the double scaling limit is a feature of 
matrix models. This enables us to extract nonperturbative information from matrix 
models. Until now, many types of matrix models have been studied in detail and 
people have made every effort to gain new information from these models. We know 
much about the matrix models which are identified to Liouville theory coupled to 
the c < 1 minimal conformal matter. We hope that we can jump over the barrier 
of c = 1 or find new universality classes in the continuum limit by modifying these 
conventional models or creating new models. 

Modified matrix models which we study in this paper have contact terms such 
as 5'(tr$^)^ in the action. The contact terms create surfaces with microscopic 
wormholes. By changing coupling constants of these square terms, we can change 
the weight for surfaces with wormholes in the matrix integral. This kind of modified 
matrix models has been studied by many people for a long time [Q]. An interesting 
point of these models is the existence of a critical value gt for the coupling constant 
g. For g < gt, the model belongs to the same universality as the g = model. For 
g > gt, as easily expected, one cannot take the double scaling limit any more and 
degenerate surfaces dominate. When we fine tune the coupling constant to gt, a new 
type of critical behavior appears in a double scaling limit and the string susceptibility 
7 jumps to a new positive value 7. At present, similar results are confirmed in the 
context of conventional matrix models — 1-matrix models, 2-matrix models and the 
c = 1 model. In general, the new exponent 7 is related to the original negative value 
7 as 



Recently, an extension of these modified matrix models is investigated [Q. In- 
stead of the square term of the puncture operator, the authors of Ref. ||^ included a 
term of the form g'{p'{0)Y to the action, where p'{0) is a representation of a scaling 
operator O by matrices. They found a critical value for g' . When we fine tune g' to 
the critical value, the gravitational dimension of this scaling operator changes to a 
negative value. 

According to our common knowledge of matrix models, the positive string sus- 
ceptibility at g = gt and the negative gravitational dimension at g' = g't are very 
strange. Recently, a simple explanation for these phenomena is proposed in lan- 
guage of Liouville theory |^ . Let us have an introductory remark for the Liouville 
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theory to explain the logic clearly. 

The path integral over a closed Riemann surface of genus g is given by the 
following form 

Zf^ = JdTj[d<i/][d(f)][db][dc\e-^^-^^-'^'-'-^''-% (1.2) 

where (p is the Liouville field and is the (p, q) minimal matter field 0. In ( |1.2D 

S"^ is the action for the matter and Sb^c the action for the ghosts. The sum of the 
Liouville action and the interaction term is written as 



S(f, + Sint 
1 
877 



Jd^z^[g^%(j)db<P - QR<P + tOmin(^) e"<^ + toO{^) e^'^]. (1.3) 



In ( |1.3| ), we have two interaction terms. The coupling constant t of the first inter- 
action term is the cosmological constant and Omin(^) is the conformal field which 
has the lowest weight /imin = Similarly, to is a coupling constant of the 

source term of the scaling operator and 0(\1/) is a conformal field whose weight is h. 
The parameters a and (3 in ( ]1.3| ) are determined by requiring the puncture operator 
P = Jd'^zO^i^e""^ and the scaling operator O = Jd'^zOe^'^ should be conformally 
invariant. The conditions for the conformal invariance are 

/imin - ^a(tt + Q) = 1, 

h-^P{(3 + Q) = 1. (1.4) 

We denote the larger solution and the smaller one by a+ and a_ respectively. Sim- 
ilarly, we have two solutions /3+ and We usually adopt a+ and /5+ as they have 
a connection to the semiclassical limit of the Liouville theory The string sus- 
ceptibility and the gravitational dimension of the scaling operator O are found by 
shifting the zero mode of the Liouville field 0, 

Q 2 
7 = 2 + -^ 



tt+ p + q - i 
d = 1-^. (1.5) 

The proposal of Klebanov and Hashimoto 0] is that the jump in the string 
susceptibility and the gravitational dimension observed in modified theories can be 



^We assume in this paper that p and q are coprime integer and p,q > 2. The central charge of 
this minimal matter is smaller than 1. 
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interpreted as a consequence of the change in the branch of a or (i. If the branch of 
a changes from (+) to (— ), the exponent jumps to a new value 



7 = 2 + ^ = —^. (1.6) 

a_ p + g + 1 



This agrees with the matrix model result This change of the branch changes 

the gravitational dimension of the scaling operator as well. The jump in the grav- 
itational dimension found in is also understood as a consequence of the change 
of the branch in /3. 



1.1 Summary of Results 

In this paper, we study matrix models modified by two square terms of the punc- 
ture operator and a scaling operator and find possible continuum limit changing the 
coupling constant of the square term of the puncture operator g and the one of the 
scaling operator g' . Our strategy to this problem is similar to 0]. We linearize the 
action by introducing baby universes and perform the matrix integration first. We 
regard the effective action derived by this matrix integration as the effective action 
of a theory of baby universes. We derive the expansion form of the effective action 
near the critical point of the unmodified theory. The free energy is expanded around 
the saddle point of the effective action and the integrations over baby universes are 
performed. We read scaling characteristics from the result of these integrations. 
Figure 1 summarizes the universality of the modified models. In this figure, the 
constants a2~^ and ^2"^ depend on the details of the matrix model. The structure 
of the universality diagram of the generic case (Fig. 1(a)) is different from the special 
case (Fig. 1(b)). As a result we find two new aspects of the modified theory. The 
first and the most important discovery is that change of the branch in /5 never occurs 
for scaling operators with gravitational dimension d> |(1 + 7). This means that we 
cannot get new critical behavior by introducing square terms of such scaling opera- 
tors. The second is that we must choose the representation of the scaling operators 
in matrices carefully to realize simultaneous tuning of the coupling constants of the 
modification terms as Fig. 1(b) i. We confirm that the tuning of the coupling con- 
stant means that the branch of the Liouville dressing a or /? changes to the opposite 
branch. In language of the theory of the baby universes, this fine tuning means 

^Klebanov and Hashimoto discussed the problem of simultaneous tuning in However, the 
fact that odd point functions of a scaling operator does not always vanish modifies their results 
slightly. 
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that the baby universe corresponding to the scahng operator becomes massless. Of 
course the relation between the modified free energy and unmodified one derived in 
is vahd if we choose the representations of scahng operators appropriately. 




Figure 1: This figure shows the universality of the modified matrix models when d < 
^(1 + 7). The universality is represented by the branches of the Liouville dressings a and 
(3 like (H — ). The left hand side is the branch for a and the right hand side the one for 
/?. When the mass matrix of the effective theory is non-diagonal, the diagram has no 

simultaneous tuning point ( ) as in Fig. 1(a). Fig. 1(b) shows that the simultaneous 

tuning is possible only when the mass matrix becomes diagonal. 

The rest of the paper is organized as follows. In section 2 we discuss the modified 
matrix models which contain the square term of the puncture operator and the one of 
a scaling operator. We find the precise form of the analytic part in the free energy 
of the unmodified model expanded around a critical point in Appendix A. Using 
this, we solve the saddle point equations and find the general form of the modified 
free energy. We read the critical exponents and complete the universality diagram 
depicted in Fig. 1. In section 3 we discuss some problems of the representation of 
scaling operators in the modified matrix models. We summarize the results and 
discuss in section 4. In Appendix A we perform matrix integration of the 1-matrix 
model with potential of the sixth order to clarify the generic form for the non-singular 
part of the free energy of the conventional matrix models. 
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2 Matrix Models Modified by Square Terms of 
Scaling Operators 



Matrix models modified by square terms of general scaling operators has studied by 
Klebanov and Hashimoto. The Ising model coupled to gravity and the generalization 
of it have been studied by them. They assumed that odd-point functions of a 
scaling operator must vanish. In general, this is not always true. This point has a 
serious influence on the universality of the modified models. Their argument must be 
modified because of their inexact assumption. Therefore we will study the modified 
matrix models carefully. 

In general the explicit form of matrix integral depends on the details of the 
model. Then we give the modified matrix models in a general way and study the 
universality of the models completely. The definition of matrix models modified by 
square terms is simply given by 

Z = jo^ e-^^ 

S = S^, + Aop{P) + E(r,)oA(0.) - ^iPiP))' - E J^iP-mr, (2.1) 

where we denote the matrix measure symbolically by D$. For 1-matrix models this 
is simply the sum over a matrix $ and for 2-matrix models it means the sum over 
two matrices. In the action, Scr is the fine tuned action of the matrix model which 
realizes Liouville theory coupled to {p, q) minimal conformal matter in the continuum 
limit i. Aq is the bare cosmo logical constant. We consider the cosmological constant 
as the coupling constant for the puncture operator P. (rj)o are the bare coupling 
constants for the scaling operators Oi which have the gravitational dimension di in 
the unmodified model. We assume that none of the scaling operators have the same 
gravitational dimension. We express the matrix representations of the puncture 
operator and the scaling operators Oi as p{P) and Pi{Oi) respectively. The coupling 
constants of the square terms g and Qi and the choice of the representations p 
and Pi characterizes the modified models. Varying these parameters we study the 
universality classes of our models in the continuum limit. We assume that the 
number of the source terms of scaling operators must be larger than the number of 
the square terms of scaling operators, i.e. , n > m. We need a source term in the 
action for each square term. 

■^We exclude the c = 1 model here. 
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From our definition of tlie modified matrix models, the jump of the string sus- 
ceptibility has been observed in the case where only the square term of the puncture 
operator exists.. In this paper, we mainly discuss the case n = m = 1. Essential 
points of our models appear in this simple case. The action of the modified model 
in this case is written as 

S = S^, + Aop(P) + rop'iO) - ^{p{P)r - J^ip'm"- (2-2) 

As the action S contains square terms we cannot perform the matrix integration 
directly. Applying a trick, we can rewrite (|2.1| ) in the following form 

Z(Ao, To, g, g') = ^^^^ jdy e"^ jdv e"^ jD^ e'^^ (2.3) 



S = S,, + (Ao - y)p{P) + (to - v)p'{0). (2.4) 
At this stage we can perform the matrix integration in (|2.3|). It is important to 



understand qualitative structure of the logarithm of the matrix integral e"^"^ 
expanded around the critical point of the unmodified model. We discuss this in 
appendix A and the result is given by 

log Jd^ e"^^ = Analytic Part + F(x,w,A^^), 

F{x,u,N^) = _ ^"^^ _ x'-^ + ■ ■ - iv°(a4 log X + ■ ■ ■) 

-A-2(a5X-2+7 + ...) + 0(iV-4), 

~N^u{bsx''+^-^ + ■■■)- iN^^ib^x^''-"' + ■■■) 

-lN\%b5X^'^-^~^ + ■■■) + ■■■ + 0{N%x'^-^). (2.5) 

It is necessary for us to know an explicit form of the analytic part of this expansion. 
The structure depends on the value of the gravitational dimension of the scaling 
operator O. On the other hand, the structure of the most singular part ( p.5|) holds 
for any scaling operators. In general the analytic part of the free energy contains 
contributions from the higher-genus surfaces. If we restrict the upper limit of the 
gravitational dimension of the scaling operator to 1, the analytic part in A and 
r comes from only planar graphs. Besides, when we limit to d < |(1 + 7), the 
non-singular part contains at most bilinear terms with respect to x and u, 

log Z{x,u) = N'^{-aix + la2x'^) + N'^{-biu + lb2u'^) - N'^boxu 

+F{x,u,N^). (2.6) 
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When we study modified matrix models, it is essential to distinguish the case of 
d < 1(1 + 7) from the case of d > |(1 + 7). As far as we consider conventional 
matrix models, the condition 3(i—l— 7 < leads < 2(i— 7 < 1 and 1 < d+l— 7 < 2. 
Thus the above form of the non-singular part is guaranteed only by the condition 
d < |(1 + 7). When d > |(1 + 7), some terms such as Ar^ appear in the bulk 
part. As far as we consider the conventional matrix models, these terms are not 
important since we discard the analytic terms when we take the double scaling limit. 
If we don't fine tune the coupling constants, these terms is irrelevant even in the 
modified theory. If we fine tune them, however, such terms cause scaling violation: 
we cannot take any double scaling limit. We will discuss this point later. For a 
while we concentrate on scaling operators which satisfy the restriction d < |(1 + 7). 

Rewriting the integration variables in terms oi x — Ao — y, u — tq — v, we find 

dx duef^^'^'l 

-00 J —00 



f{x,u) = -—[Ao^-2x{Ao-aig)+x^{l-a2g)]-NXux 
^9 



'^9 



-[To' + 2u{ro - hg') + u\l - hg')] + F{x, u, N^). (2.7) 



We keep the leading non-singular parts in f{x, u). This analytic part consists of the 
terms which are bilinear at most with respect to x and u. These are necessary for 
applying the expansion around the saddle point explained later. 

We can regard the function f{x,u) as an effective action for a theory of baby 
universes. The fields x and u effectively correspond to the operators that make a 
puncture to a surface. This puncture is characterized by the scaling operator. Thus 
we distinguish the baby universe of the puncture operator from the baby universe 
of the scaling operator O. The bilinear terms of f{x, u) are related to the masses 
of these baby universes. We have two baby universes associated with the puncture 
operator P and the scaling operator O. However the identification of these two baby 
universes is a little complicated. To make pure baby universes, it is convenient to 
rewrite /(x, u) in the following form. 



f{x,u) = AT^ 



Al Tl 

9 9' 



X 




X 


T r 


u 




u 





+F{x,u,N'), 



Oo 

9 ^ 



bo 



9' 



bo 




X 


-b2 




u 



(2.8) 
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where Ai = Aq — aig, ti = tq — big'. We drop irrelevant constant terms from the 
effective action. The mass matrix of the baby universes 



- a2 



bo 



[2.9) 



is diagonal only when = 0. Then for bo 0, x and u are not eigenfunctions of the 
mass matrix. When bo = 0, the non-singular part has already diagonalized and we 
can regard x and u as baby universe fields. Whether 60 = or not depends on the 
choice of the representations p and p' in ( |2.2| ). This parameter bo is nothing but the 
one which distinguishes the generic case (Fig. 1(a)) from the special case (Fig. 1(b)). 
We must discuss these two cases separately. 



2.1 60 7^ 



For 60 7^ we diagonalize the non-singular part by redefining w = u + ^z^^^;. This 
redefinition is valid for g' ^ 62^^. For a while we assume g' 7^ 62"^ and consider the 



exceptional case g' = b2 ^ later. Redefining the coupling constants as 



(Ao - aig) 

hg', 



bo9 



To 



1 - b2g 



7(7-0 - big'), 



(2.10) 



we find 
f{x,u) 



f{x,w 

iV2 



A T_ 

9 9' 



X 




X 


T " 


w 




w 





a2 - 



1-629' 



+F{x,w,N^). 



h'b2 





X 




w 



(2.11) 



The leading singular terms in F{x, w, N"^) are the same as the ones in F{x, u, N"^). 
The redefinition of the baby universe for the scaling operator O keeps the leading 
singular terms of the effective action. When we seek for the double scaling limit of 
the modified model, we must take care of the form of the singular part. When we 
diagonalize the mass matrix of the baby universes, the large N expansion form of the 
singular part must be written as ( p.5|) characterized by 7 and d. This is the reason 
why we redefine only u. The mass matrix can be diagonalized by redefining x. 
However, this changes the form of the leading singular part. 



9 



As the non-singular part is diagonal in x and w, we can read the masses of these 
baby universes from each diagonal component, 



, --a2-^. (2.12) 

From assumption w baby universe cannot be massless. The last contribution to the 
mass square of x baby universe comes from the two point function of the puncture 
operator and the scaling operator O. We perform a large expansion around the 
saddle point {xs,Ws) which is given by the equation = |£ = 0. The location of 
the saddle point is closely related to the masses of two baby universes. Thus we 
discuss separately by the sign of the mass eigenvalues 

First we consider the case g' > 62"^- This means that w baby universe is tachy- 
onic. This case corresponds to the region above the dotted line in Fig. 1(a). It is 
easy to see there is no stable saddle point near {x,w) = (0,0). The large expan- 
sion of the free energy which is found by expanding /(x, w) around the saddle point 
cannot be written in the form of ( |2.5| ). Thus this case causes scaling violation and 
we cannot take any double scaling limit. We can, however, take a continuum limit 
in which only planar graphs survive. 

Next we consider the case g' < &2~^ in which w is massive. The region below 
the dotted line in Fig. 1(a) corresponds to this case. We divide this case further into 
three: the x baby universe is massive, tachyonic and massless. 

First the tachyonic case nix^ < causes scaling violation for the same reason as 
the tachyonic w case. In general when at least one of the baby universes is tachyonic, 
we cannot take a double scaling limit. 

The case in which both of the baby universe are massive corresponds to the 
inside region of the hyperbolic curve in Fig. 1(a). The solution of the saddle point 
equation in this region is 

A 



Xq 



0{A 



l-7\ 



Shifting the integration variables, x 
log Z {A, T,g,g') 



x — 



r 

g'mj 



w 



0(A' 



W — 



(2.13) 



we have 



log / Ndx 



X exp 



F\x + 



Ndw exp 

o 

A 



X 



w 



w + 



T 



g'm^ 



(2.14) 
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We rescale integration variables as x = xN"^/^"^ '^^a^'^^^'^ w = wN"^^^ "^^Z*^^ '^^ and 
choose t = ^^A^2/(2-7) 1/(2-7) to = ^^iV2(i-'^)/(2-7) as renormalized scaling 

variables. Then we have 



logZ(A,T,5(,5(') 

POO 

= log / A^-^/(2-7)^£exp 



-^-27/(2-7). -2/(2-7)^ 2^ 



^— oo 

X exp[F(a; + t,w + to)]- 



^^2(2d-7)/{2-7)^^2^2 



^;2 



(2.15) 



We keep t and to finite and take the large limit. Then the integrand becomes a 
product of delta-functions S{x)6{w). We have finally 



\ogZ{A,T,g,g') = F{t,t 



(2.16) 



This scaling function F{t, to) is the same as ( |A.14D appeared in the g = g' = 
theory. Thus the modified theory belongs to the same universality class of the 
unmodified theory. 

In the massless case rria;^ = 0, the location of the saddle point has a new scaling 
nature, 



Xq 



\ +0{A—), Ws 



90-3 



g'm.. 



0{A'^+^-^). 



(2.17) 



We shift one of the integration variables, w = w r^-^. The gaussian term of w 

becomes a delta- function 5{w) in the large limit. Thus we can simply replace w 
in the integrand by 2 . After this replacement we have 

/oo 
Ndxef^^\ (2.18) 
-00 

N^A 



fix) 



-X + F \ X, 



g'myj'^ 



(2.19) 



The problem reduces to a theory of a single baby universe. We expand the effective 
action f{x) around the saddle point Xs and find the A^^-order contribution 



log Z {A, T,g,g') 



2-7 



«3 



1 /(1-7)AV- 



2-7 



1-7 



ga3 



+ • 
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d+1-7 
1-7 



+ 



1 - 7)A\ 



+ 



(2.20) 



This form of the large N expansion is similar to ( |2.5| ) if we take new string suscep- 
tibility 7 and new gravitational dimension d of the scaling operator as 

7 = ^. (2.21) 

7 — 1 7 ~ 1 

These new values can be understood as the result of the change of the branch in 
a in the Liouville interaction e"*^ from to a_. For the string susceptibility, we 
explained it in §|T]. For the gravitational dimension, we can easily check d = \ — ^ 
agrees with the above result. The branch for (3 is unchanged. 

The massless x region lies on the hyperbolic curve in Fig. 1(a). We consider this 
curve as a critical line. When we fix the value of g' , the location of the critical point 
Qt varies with g' . Of course, if we turn off the square term of the scaling operator 

0, gt = a2-\ 

It is easy to see that the double scaling limit is realized by choosing new scaling 
variables as 

i = ^iV2/(2-7)((l_^)a3)-(l-^)/(2-7)^ 

g 

tQ = ^ AT2{l-d)l(2~^) ^ ^ ^Y2(l-d)/(2-7)_ /2.22) 

In the double scaling limit the free energy acquires new scaling nature 

\ogZ{^,r,g,g')=F{t,to), (2.23) 

+2—7 oo oo oo _ 1. 

ni to) = 7— -«i logt-E -aHt^'-'^^'-'^-i: E bnM'^'-'^^'-'^ {tot'-') . (2.24) 

^ h=2 h=Ok=l 

We can derive a relation between the modified free energy F and the unmodified 
one F appeared in ([A.14|) . To do this we rescale the integration variable in ( p. 18 



as 



t = xN"^/^"^ '''^as^/*^^ and take the double scaling limit. We find a formula similar 
to the one derived in 

/oo _ 
e^**+^(*'*'^\ (2.25) 
-oo 
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where e is a constant which depends on only 7 (or 7), e = (1 — 7)'-^ ^-'/^^ . This 
relation can be considered as a nonperturbative definition of the new scaling function 
F. 

Now we come back to the case g' = 62"^ we have left at the beginning of this 
subsection. The baby universe w is not well-defined for this case. As there is 
no quadratic term of u in the effective action f{x,u), we cannot absorb the term 
—N'^boxu by redefining the baby universe u. For finite g, the result of the large 
expansion around the saddle point cannot be written in the form of (|2.5| ) by any 
means. This means that we cannot take a double scaling limit in this case. We 
cannot regard the modified theory in this parameter region as an effective theory of 
baby universes. When g = which means the baby universe corresponding to the 
puncture operator is absent or extremely heavy, we can find a new scaling nature in 
the theory. To show this we return to the starting point, 

/oo 
due^^^\ (2.26) 
-00 

f{u) = N%u (^0 - I - |ao^ + F(Ao, u, N'). (2.27) 

As we have no quadratic terms of u, the u baby universe is massless. Defining new 
coupling constants as A = Aq, r = tq — |^ — f^Ao, we find the location of the 
saddle point 

Us = -^rA-2'^+-^ + 0(A^-''). (2.28) 
64 

This estimate of the saddle point is valid for sufficiently small A. Integrating around 
Us and choosing the scaling variables as 

t = AAr2/(2-7)«3i/(2-7)^ to = &2rA2(i-'^")/(2-7)^ J = 7 - rf, (2.29) 
we find the modified free energy in the double scaling limit 

\ogZ{A,r,0,b^') = F{t,to), (2.30) 

J.2— 7 00 

Fit,io) = -j-o r-ailogt-Ea/.^^'-"^^'-"^ 

(2-7)(l-7) f^^ 

00 00 _ {, 

- E E V.^^'^^^^'^'H^ot'^'-O ■ (2-31) 

h=Ok=l 

Thus the string susceptibility of the modified theory is the same as the one before 
modification. On the other hand, the scaling operator O acquires a new gravitational 
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dimension d = -y — d. In language of Liouville theory, this is explained by the change 
of the branch in the Liouville dressing e^"^, that is, 

d = j-d=l~—. (2.32) 

This universality appears at a single point {g,g') = (0,62^^) in Fig 1(a). In a sense 
we can say the value g' = 62"^ is special. Thus we may call this special value of g' 
as a critical point for g'. 

Rescaling the integration variable as Iq = 'uA^2(i-c()/(2-7)^ ^j^g^^ ^j-^g mod- 

ified free energy is related to the unmodified one 

/oo _ 
dto e*«*«+^(*'*°). (2.33) 
-00 

This relation resembles ( p.25| ). Now the integration is performed over to not t. 

So far we have studied the generic case 7^ 0. It is worth noting that for bo ^ 0, 
the change of the branch in a and (3 never occurs at the same time. This is deeply 
related to the fact that we cannot make the two baby universes corresponding to 
the puncture operator and the scaling operator massless at the same time. 



2.2 60 = 

Now we study the special case bo = when the mass matrix of baby universes is 
diagonal. The universality diagram for this special case is shown in Fig. 1(b). A 
remarkable difference from the generic case is that it is possible to fine tune g and g' 
simultaneously. This means that the modified theory may have another new scaling 
nature which is present in the generic case. The simultaneous tuning corresponds 
to the single point {g,g') = (o.2~^i^2~^) in Fig- 1(b). The work of Klebanov and 
Hashimoto covers this special case. The non-singular part of the effective action 
f{x,u) of baby universes ( |2.8| ) has already diagonalized. We can treat x and u as 
pure baby universes. We interpret that x corresponds to the puncture operator and 
u to the scaling operator O. The mass squares of these baby universes are read from 
the effective action, 

1 1 

m^'^ = 02, mj = — - 62. (2.34) 

9 9 

As in the generic case, we divide into five cases by the sign of the mass squares. In 
most cases the analysis is common to the generic one and the answer is the same as 
IQ. We do not repeat in detail and show only some essential results. 
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First we consider the case where x and u are massive. This case corresponds to 
the inside region of the rectangle in Fig. 1(b). In this case the integral becomes a 
product of delta-functions around the saddle point in the large limit. This allows 
one to replace x with Xs = j^;^ and u with Us = jZb^i where A = Aq — aig, 
T = To — big'. Therefore the modified theory is in the same universality class as the 
unmodified one. We select the renormalized scaling variables, 

^ -iV2/(2-7)a3V{2-7)^ t^ = _L_N'^i-m2-^)_ (2.35) 



1- a2g 1 - b2g' 

In the double scaling limit, the sum over surfaces converges to F{t,to), the gener- 
ating function of the g = g' = theory. 

The second case is 5^ = 02"^, < g' < ^2"^, ^-e. , the case of massless x and 
massive u. Contrary to the generic case, the massless condition nix'^ = requires 
that the coupling constant g be fine tuned to 02"^. The double scaling limit is 
realized by choosing scaling variables as 

A = Ao - — , T = To - big', (2.36) 

^2 

t = a2A2/(2-7)(i _ ^)a3-(^"^)/(2-^), to = ^^iv2(i~d-)/(2-^)^ (2.37) 

1 - b2g' 

where 7 and d are the new string susceptibility and the new gravitational dimension 
of the modified theory, 

7=^. rf=^. (2.38) 

7 — 1 7^1 

These values are explained by the change of branch in a. The massless condition 
for X baby universe is nothing but the fine tuning of the coupling constant g to 
its critical value gt = a2~^- As the consequence of this fine tuning, the branch 
in a changes. In the double scaling limit, the free energy becomes F{t, to) which 
appeared in ( p.24|) in the generic case. Thus the modified theory belongs to the 
same universality class as the case m^.^ = and m^^ > in the generic case. 

The third case is < g < 02^^ and g' = b2^^,i-e. , the case of massive x and 
massless u. When 60 7^ 0, the coupling constant g must be zero to take a double 
scaling limit. In this case, however, g is not restricted to zero. The scaling variables 
are defined by 

A = Ao - aig, r = tq - ^, (2.39) 

02 

-iV2/(2-7) ^^1/(2-7)^ ^-^^^^^^2(l-d-)/(2-7), (2.40) 



0-29 
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The scaling operator O gets a new gravitational dimension d = j — d. This new 
gravitational dimension is understood by the change of the branch in jS. The scaling 
function F{t, O) is the same as ( |2.31| ) in the generic case. 

The fourth case corresponds to simultaneous tuning of the two coupling constants 
g and g'. The solution of the saddle point equation is 



'1- 



Xc 



03 



2+7 



(2.41) 



where A and r are given by 



A = An- 



02 



r = To 



61 



(2.42) 



If we denote the derivative of the function /(x, u) with x by /'(x, u) and the deriva- 
tive with u as /(x, u), we have 



-1 u -1^ 



logZ(A,r, a2-\h 



f"f - if ) 



+ 0{N 



-2\ 



(2.43) 



-I X=Xs,U=Us 



After some calculation we find new scaling nature in the free energy. The double 
scaling limit is realized by selecting the scaling variables as 

i= a2A2/(2-7)(i _ ^)a3-(i-^)/(2-7)^ ^ ^^^^2(l-d-)/(2-^)^ ^2.44) 
The new string susceptibility 7 and the new gravitational dimension d are 

7 



7 Q 
= 2 + — 

7 — 1 Q!_ 



d 



d 



7-1 



1 



(2.45) 



These new values indeed correspond to the simultaneous change of the branches in 
a and (3. Keeping these two scaling variables finite and letting to infinity, we find 
the scaling function of the form, 



\ogZ{A,T,a2-\b2-')=F{t,io), 



(2.46) 



+2—7 00 00 00 _ u 

F{i, to) = ^--ai logt-E -ani^'-'^^^'-'^-T. E hj^'-'^^'-'^ (tot'-') . (2.47) 

h=2 h=0 k=l 



2-7 
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The scaling function ( p.47| ) may be expressed by an integral over scaled coupling 
constants, 

t = xiV2/(2-7)a3l/(2-7)^ to = uN^^^""^ ' (2.48) 

We find a simple relation 

/oo roo _ _ 

dt / dto e'^*+^oto+F{t,to) ^ (2.49) 
-oo J —oo 

where e = (1 — ^)(i~t')/(2-7) jg ^ constant which depends on the universality of the 
unmodified theory. The appearance of this new universality class is a feature of the 
60 = case. 

The last case is the one where at least one of the baby universes is tachyonic. 
This corresponds to the outside region of the rectangle in the universality diagram 
Fig. 1(b). As there are no stable saddle points near x = m = 0, we cannot take any 
double scaling limit. Thus we are not interested in this case. 

2.3 Scaling operators with d > |(1 + 7) 

Finally, we briefly comment on the case d > ^(1+7). In this case, the effective action 
of baby universes cannot be written in the form ( |2.8| ). There are some extra higher- 
order terms in the non-singular part. Nevertheless when {1 ~ a2g){l — h2g') — bo^ g g' > 
0, we can take a double scaling limit by choosing the scaling variables as 

t oc AN'/^'-'\ A = Ai--^n, 

1 - b2g' 

to oc riV2(i-d)/(2-,) ^ ^ = ^^__^Ai, (2.50) 

1 - «25' 

where Ai and ti are appeared in ( p.8|) . In the double scaling limit, the modified 
theory is in the same universality class as the unmodified theory. Such a double 
scaling limit is also possible when d < |(1 + 7). The condition (1 — a2g){l — 
b2g') ~ bo^gg' > corresponds to the case both x and w baby universes are massive. 
Therefore we have two different ways of taking the double scaling limit when (1 — 
a2g){l — b2g') — b^^gg' > and d < |(1 + 7). We cannot take a double scaling limit 
out of this region unless we turn off the coupling constant of scaling operators with 
(i>|(l + 7). 

The reason why we are not interested in the scaling operators whose gravitational 
dimension satisfies d > |(1 + 7) is that we cannot change the branch in o; or /3 in 
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a double scaling limit. The square term of such operator is irrelevant when the 
baby universe associated with the operator is massive. Otherwise it destroys scaling 
property of the theory in the continuum limit. 

3 Representation of Scaling Operators and the 
Modified Matrix Models 

In appendix A we have studied the 1-matrix model with potential of sixth order. We 
have learned that both tr and tr + q;$^), a ^ —15 behave as the puncture 
operator in the continuum limit. This example teaches us that the way to represent 
scaling operators by matrices is not unique. For a while we discuss the problems of 
the representation of the scaling operators. 

In the 1-matrix model of appendix A, a special linear combination tr — 15$^) 
behaves as the scaling operator with dimension |. In general, we can produce 
an unknown matrix representation of a scaling operator like this way. Prom two 
representations pi , p2 of a scaling operator with gravitational dimension d, we can 
make a new representation p = pi + rp2 so that it represents a scaling operator 
with dimension larger than d by choosing r appropriately. It is essential that the 
new representation is always a linear combination of the two representations. This 
fact comes from the universality of the matrix models, that is, the form of the 
scaling function is unchanged up to the dilatation of the scaling variables. Thus 
we can cancel the leading singular terms in the one point function of these scaling 
operators (pi), i = 1, 2 by defining a linear combination p — pi + rp2- 

In the last section, we have learned that the scaling nature of the modified matrix 
model depends on not only g and g' but also bo. The value of bo is determined 
by the choice of the representations p{P) and p'{0). Is it be possible to choose 
these representations so that bo vanishes? The answer of this question is yes. The 
two point function {p{P)p'{0)) of the unmodified model determines the value of 
bo- Even if bo 7^ 0, we can find a new representation of the scaling operator O 
that makes the mass matrix of the baby universes diagonal. To do this, we need 
another representation of the scaling operator p"{0) such that the two point function 
{p{P)p"{0)) does not vanish. Such a representation always exists. We define a new 
representation of the scaling operator by a linear combination of these, p — p' + rp". 
If we replace p' by this p, the parameter bo must be a linear function of r. Thus we 
can make 60 = by appropriately selecting r. This means that we can always fine 
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tune g and g' simultaneously for the scaling operator whose gravitational dimension 
satisfies the condition d < |(1 + 7). 

4 Discussion 

In this paper we have confirmed that the string susceptibility and the gravitational 
dimension in the double scaling limit are understood by the idea that a fine tuning 
of square terms changes the branch of Liouville dressing of some scaling operators. 
This fine tuning is possible when the mass matrix of the baby universes is already 
diagonal. Even if it is not diagonal, we can diagonalize it by redefining the baby 
universe corresponding to a scaling operator except the puncture operator. When 
the mass matrix is diagonalized by this procedure, the simultaneous tuning of square 
terms is impossible. Our new discovery is that the representations p{P) and p'{0) 
have a effect not only in the critical values for g and g' but also in the possibility of 
simultaneous change of the branch in a and f3. 

The existence of an upper bound in the gravitational dimension of the scaling 
operator above which we cannot fine tune the coupling constants of square terms is 
interesting for us. This upper bound means that the branch of the Liouville dressing 
for such operators cannot be changed in the modified theories we have studied. If 
we fix the type of the minimal conformal matter coupled to gravity, the number of 
allowed operators is always smaller than the number of the gravitational primary 
fields. We think this is related to the fact that the difference between the behavior of 
g/3+<^ and e^~'^ becomes large for such operators in the weakly interacting region of 
large in Liouville theory. Our result doesn't insist that it is impossible to change 
the branch of such operators by any means. It may be realized by another type of 
modified matrix models. 

We have restricted our argument to the case where there are at most two square 
terms, i.e. , m = 1 in ( |2.1| ). Of course we can generalize to the case m > 2. 
Then the size of the mass matrix of the baby universes becomes (m + 1). The 
simultaneous tuning of (m + 1) square terms is possible only when the mass matrix 
is diagonal and all the scaling operators satisfy the condition d < |(1 + 7). Choosing 
the representations of scaling operators appropriately, we can always make the mass 
matrix diagonal. The idea that a fine tuning changes the branch of Liouville dressing 
of some scaling operators is valid as well. It is interesting that the matrix models 
we have considered and the theory of baby universes are intimately related like this. 
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When the square terms are fine tuned, the continuum theory doesn't have semi- 
classical limit. This fact makes the modified matrix models more complex. However, 

the results of the modified theories strongly suggest the existence of such a strange 
Liouville theory. We hope that physical aspect of these fine tuned theories becomes 
clear. Our work will be easily applied to the c = 1 matrix model. 
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Appendix A General Form of the Results of the 
Matrix Integration 

To study the modified matrix models, it is necessary to know a general form of 
the free energy expanded around a given critical point. In this appendix we study 
a pedagogical example and deduce the general form of the free energy from this 
example. We use the 1-matrix model with potential of the sixth order as the example. 
This is the simplest model if we want to include a source term of a scaling operator 
with gravitational dimension greater than 0. The partition function is defined by 
the integration over an N x N hermitian matrix $, 

where A and k are parameters. Usually we relate A with the cosmological constant in 
the Liouville theory. It is well-known that the large- AT expansion of the free energy 
corresponds to the topological expansion. 

Using the method of orthogonal polynomials or the steepest descent method, we 
can find the contribution of planar graphs, 

(A, k) = -36A«;a^° + {15k- 6X'^)a^ + (3A + 20K)a^ + (^5A - ^) 

-^«' + ^ + ^loga^ (A.3) 
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where a satisfies an equation 



GOna^ + 12\a^ - + 1 = 0. (A.4) 
From ( [A .31 ) and ( [A .41) we eliminate A and find 

e(°)(A, k) = ^(1 - a^){9 - a^) + ^ loga^ - 7Ka^ + 2Ka^ + 30/c^a^^. (A.5) 

This free energy shows a singular behavior if we fine tune A and k. The critical 
point such that ( |A.4| ) has triple root in is known as the = 3 critical point of 
1-matrix model. At this critical point, the theory is identified with the Liouville 
theory coupled with (2,5) minimal conformal matter. The /c = 3 critical point is 
achieved by fine tuning the parameters to the values. 



_ 1 , 1 

'^'""1620' 36' 



A^ = — , a? = 3. 



To find the singular behavior for the free energy near the critical point, it is conve- 
nient to introduce = a^(l — 6), A = Ac — A. A is nothing but the cosmological 
constant. We need the expansion of 6 by As up to the seventh order, 

b = 3(4A)3 - 6(4A)3 + 9(4A) - 10(4A)3 + 7(4A)3 (4A)3 + ^(Aa). 

3 

Defining r = — k, we expand e*^°^(A,r) around A = r = 0. The r-independent 
contribution is 

e(°) A,0) = + -log3 + 0(b^) 

^ ' ' V 40 2 ^ y 20 10 20 40 14 ^ ' 

19 1 , „\ 27 . 2 19683 .... 7 



-40 + 2^°s2j-T^ + ^^^^^ 
+C(At). (A.6) 

As the leading singular term behaves as A^""^, the string susceptibility is — |. This 
exponent suggests that the continuum limit of this matrix model is identified with 
the Liouville theory coupled to (2, 5) minimal conformal matter. 

The singular behavior for the linear and quadratic terms in r is derived by 
differentiating with respect to r 



= -27 + 34992A - 98415(4A)3 + 0(A3), (A.7) 

T=0 
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52g(0) 



= 437400 - 7873200(4A)3 + 0(A3). (A.8) 

r=0 

Gathering these contributions, we arrive at the final form of the free energy, 



27 1 1 

( A, r) = A + - • 2196A2 - 27t + - • 437400r2 + 34992At 

5 2 2 

19683 (^A)i + ----98415(4A)tr + --- 



28 

•7873200(4A)5t2 + --- 
+0(A-ir3). (A.9) 

From the definition of the gravitational scaling dimension, one can easily find that 
tr corresponds to the puncture operator [d — 0) in the continuum limit. 

It is well-known that in the (2,5) minimal matter theory there are two primaries: 
the identity operator and a operator with conformal dimension — |. After coupling 
with gravity, the latter becomes the puncture operator and the dressed identity 
operator gain the gravitational dimension |. We want to replace tr $® by the dressed 
identity operator. The operator with d = \ should be represented by a linear 
combination of these operators. Thus we replace tr by a trial operator tr — 
tr + q;$^) and determine a so that behaves as the scaling operator with 
d—\. To this end we consider a matrix integral. 

It is not necessary to repeat calculation. What we have to do is just replacing A with 
A + a(^K — Kc) in the previous result. As the location of the critical point doesn't 
move, the definition of the scaling parameters A and that of r are the same as 
before. After some calculations we find the expansion form of e°(A, r). 



=(0) 



(A,r) 

: -^A + ^ • 2196A2 - 27 + l) r + ^ ■ 2196(q;^ + 24q; + 150)t^ 

+2196(q; + 12)Ar - — (108)i A^ + • • • 
28 



+r 
1 



^(a + 15)(108)i At + i f^a + 48) (108)iAl + • • • 



3' ' 2 V 5 ; 

2 + 15)^(108)3 A3 + ^(a + 15)(13a + 285)(108)i As 
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-^(108)^(2a2 + 90a + 945)A + ■ ■ 

+0(A-ir=^). (A.ll) 

If we set a = —15, the terms proportional to A^r, A^r^ and A^r^ vanish and the 
operator tr (9_i5 = tr ($^ — 15$^) behaves as the operator desired. It is worth noting 
that this e^°-'(A,r) contains odd-order terms in r which is absent in the paper [Q. 
They assumed that odd-point function of the scaling operator vanish. Our example 
shows that this is not generally the case. 



Our result ( |A.ll ) offers a general form for the generating function of the corre- 



lation function of a scaling operator O. We consider A as the cosmological constant 
and r as a coupling constant of the scaling operator and define the matrix integral 
Z{A,t) as (|A.10|) . We fine tune parameters of the matrix model so that the string 
susceptibility becomes 7. We choose the matrix representation of O so that the 
gravitational dimension becomes d at this criticality. In general the free energy can 
be written in the following form 

log Z(A, r) = Analytic Part + F{A, t, N^), 

FiA,T,N') = N'(^-j^—^^A'~'' + ..^-N%a,logA + ...) 

-N-\a,A-^+'' + ---) + 0{N'^), 
-N^T{h-^A'^+^-< + ••■)- \N\^{hiA^'^-^ + ■■■) 
-iN^T^hA^"^-^-^ + ■■•) + ■■■ + C(A^VA^-I). (A.12) 

Here we separate non-singular part from the singular part F{A,t, N"^). The struc- 
ture of the singular part is indeed consistent with the result our example ( |A.11| ). 



This scaling property of the singular part comes from the fact that one can take a 
double scaling limit. The double scaling limit is achieved by choosing the renormal- 
ized scaling variables as 

t = AAr2/(2-7)^3l/(2-7)^ ^ ^^2(l-d)/(2-7). (A.13) 

In the double scaling limit with fixing these scaling variables finite, the singular part 
of the free energy becomes 

j-2— 7 00 

(2-7)(l-7) f^^ 

k 



h=Ok=l 
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This scaling function depends only on the scaling variables t and to and is the sum 
over the contributions of the surfaces with different topologies. When we take the 
double scaling limit, we discard the analytic part. Therefore its explicit form is not 
necessary in the conventional matrix models. However, when we modify the matrix 
models by square terms, this non-singular part becomes very important as discussed 
in §g 
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